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1 Introduction 

The term PT-symmetric quantum mechanics, although defined to be of a much 
broader use, was coined in tight connection with C. Bender's analysis of one- 
dimensional Schrodinger Hamiltonians with potentials 

V{x)^x^{ixy, (1) 

see [1] (we will call them Bender Hamiltonians in the following). This class 
of operators - characterised by one real parameter e - contained the probably 
most exploited Hamiltonian in the history of physics: the e = harmonic os- 
cillator; but on the other hand, the other members of the family were strange 
Hamiltonians with imaginary potentials which do not appear physical at all. 
The aim of the suggested 'PT-symmetric treatment was to give them an accept- 
able interpretation and it became a rather standardised procedure since then. 
One of the key facts which allowed for the physical interpretation was that the 
spectrum of these operators were real at least for some non-zero values of e, a 
fact which stems from observation that first, for 'PT-symmetric operators the 
complex eigenvalues appear in mutually conjugated pairs, and second, that the 
spectrum depends continuously on e. Thus the boundary between real-energy 
and complex-energy domains can lie only at an exceptional point, i.e. a point 
where at least two eigenvalues merge. Consequently one has to depart at least 
some finite distance on the e-axis from the non-degenerate Harmonic oscillator 
to see any selected eigenvalue complexify. What actually happens in the dis- 
cussed case is visible at the classical figure Fig. [TJ the spectrum is real for all 
positive e, while at negative e the lower the energy is the longer it remains real. 

However, the validity of the continuous energy dependence assumption is not 
obvious. An immediate question of the reader could be: What happens when 
e — 2, which yields an obviously unphysical —x'^ potential with unbounded 
continuous spectrum? It is well known now, that the correct response has to 
deal with boundary condition in the complex plane. The Hamiltonian which is 
the proper continuation of the Harmonic oscillator in the parameter e is defined 
on space of functions which are square-integrable not on the real line, but on 
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Figure 1: Well-known dependence of the eigenvalues of the Bender Hamiltonian 
on e, non-tobogganic case. Complex energies are not shown. Based on numerical 
calculation it is expected that the value at which the n-th energy becomes 
complex, tends to zero as n — > oo. 

some asymptotically straight contour in the complex plane which lies in the 
correct Stokes sector of the Schrodinger equation. Here, correct means that 
the sectors are itself a continuation of those sectors which contain the real line 
in the £ = case. These sectors (also called "wedges") are turning down in 
the complex plane as e increases and for s > 2 they no more contain the real 
line. Therefore, the conventional — Hamiltonian on L2{M.) is not the only 
Hamiltonian which deserves this name; in some sense its analogue defined by 
the same differential equation but with complex boundary conditions is more 
natural. 

For sake of clarity let us recall that for potentials we are dealing with the 
exact choice of the integration contour is irrelevant as long as it lies inside 
the Stokes sector asymptotically. It became customary to omit mentioning the 
contour altogether, speaking only about boundary conditions imposed inside a 
particular sector. On the other hand, the contour is a convenient means for 
defining the scalar product and, practically, some concrete choice of the contour 
is necessary for numerical computation; thus we will speak about integration 
contours rather than wedge-defined boundary conditions in the rest of the paper, 
keeping in mind that due to the potential's analyticity distinct contours with 
same- wedge asyniptotics are equivalent It may be also noted that in the special 
case £ = 2 it is enough to pose the boundary conditions on the shifted real line, 
i.e. X — t — ic with arbitrary positive constant c (i is a real parameter). For 

^In fact one can even release the asymptotic straightness condition provided that the 
contour does not oscillate too rapidly in the asymptotic region. The technical details of 
equivalence between contour integrability and boundary conditions in infinity arc beyond the 
scope of this article. 
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higher e one has to use bent contours to remain inside the Stokes wedge. 

The existence of aforementioned continuation is interesting from different 
points of view. For example, it is relevant to the Dyson argument about con- 
vergence of the perturbative series, which is roughly stated as follows: having a 
potential + gx"^ (or its field-theoretical counterpart, the argument was orig- 
inally formulated for quantum electrodynamics [5]) one may use perturbative 
expansion in 17, but since for any negative g the spectrum collapses to (—00, 00) 
and the perturbative calculation cannot produce this, the convergence radius of 
the series ought to be zero and the series is at best asymptotic. The existence 
of a negative-g Hamiltonian with discrete real and below bounded spectrum 
invalidates the argument's core assumption: it is now entirely possible that the 
series converges and gives the spectrum of the complex-plane —x"^ potential 
when evaluated at negative g, viz. [3]. 

2 Quantum Toboggans 

The observation which we want to point out now is that to force the integration 
contour to be asymptotically straight and inside the correct Stokes wedges is 
insufhcient to uniquely define the spectra of Hamiltonians with ([!]) . One has to 
take into account that for non-integer e the potential lives on multiple Riemann 
sheets and it matters how the integration contour is distributed upon the sheets. 
In our presently discussed class of Hamiltonians, the only singularity lies at 
a: = 0, which allows us to use only one winding number A to fully characterise 
the contour. Hence instead of one PT-symmteric continuation of the harmonic 
oscillator we obtain an infinite series numbered by distinct integer values of 
A. The A = trajectory represents the only usually considered case. The 
higher-A contours define distinct Hamiltonians which are sometimes referred 
to as quantum toboggans, the reason of such denomination is clear when one 
imagines the Riemann sheets forming a helicoicH. 

A technical note: To take profit from the fact that real eigenvalues have to 
merge before becoming complex, one must take care about the PT-symmetry. 
In particular it dictates the position of the branch cut. It has to be chosen to lead 
upwards along the positive imaginary axis. On the other hand, the potential 
has to be PT-symmetric on the contour, which, for A 0, forces the contour 
to pass below the singularity. If one wanted to have the contour bypassing zero 
from above, which is an equally reasonable option for the harmonic oscillator 
where no singularity is present, one has to change the potential accordingly 
to preserve the PT-symmetry of the potential on the contour. Without much 
surprise this yields 

V = x^{-ixY (2) 

Such change of potential together with the change of the contour would obvi- 
ously be nothing than the reparametrisation such it leaves the spec- 

^ Suggested imagination is, of course, consistent only when infinite nuber of sheets are 
present. 
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Figure 2: Considered integration contours: a straight one (winding number A = 
0) is drawn in the dashed hne, while its most elementary tobogganic counterpart 
(A = 1) is depicted in solid line. The latter contour is consistent with VT- 
symmetry only if the the branch cut aims downwards, therefore we need to use 
an inverted one with potential ([T]). 
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Figure 3: "PT-symmetric contours for A = 1,2,3 of the potential x^iixy. In 
contrast with Fig. [2] the A = 1 contour is chosen inverted with respect to the 
real axis, conforming the PT-symmetry assumption. 

trum intact. We can similarly disregard the differences within analogous mu- 
tually conjugated pairs of trajectories and potentials even in more complicated 
settings. Therefore, without loss of generality, we will stick to the definition ([T]) 
and keep in mind that also the contour must conform to the PT-symmetry. 

The integration contours are schematically drawn in Fig. [2] and Fig. [2l It 
turns out that the choice of the A = 1 contour in Fig. [5] is inconsistent with 
PT-symmetry, instead one has to use the up-down inverted curve, as depicted 
in Fig. [2] A general rule is that the center point of the contour lies below zero 
(for concreteness say sX x = — i) on the principal sheet, which is the one where 
y(— i) is real. 
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3 Numerical Analysis 



Unfortunately, the considered potentials are not exactly solvable. Therefore we 
have to rely solely on numerical treatment. We use a relatively straightforward 
method for computation of eigenvalues; we are interested whether real eigenval- 
ues are present - the existence of complex part of the spectrum can be partially 
deduced from the distribution of the exceptional points - this allows us to signif- 
icanly simplify the calculation. First of all, a starting energy is chosen. Then we 
calculate two independent solutions of the differential equation with the chosen 
energy; let us denote them ipi and 'ip2- The equation has been solved on the 



contour parametrised as 

X = t — i for A = 0, (3a) 

X = ie'*e(7r - |t|) + (< - TT - i)e(|i| - tt) for A = 1, (3b) 

X = -ie'*e(27r - |i|) + (t - 27r - i)e(|t| - 27r) for A = 2, (3c) 



and so on; i.e. the contour consists of a circle of unit radius A-times encircling 
the singularity in a- = and a straight line parallel with the real axi^ matched 
to the circle at a; = — i. For concreteness the initial conditions for ipi^2 are set 
in the centre of the contour (i.e. t = 0) to satisfy 

V'i(0) = 0, 7/^i(0) = l, (4a) 
^^2(0) - 1, i^'M - 0. (4b) 

The number E is an eigenvalue if there exists a linear combination of ipi and 
'ip2 which is integrable. Because the asymptotics of the solution is exponential, 
this is equivalent with the existence of a linear combination tending to zero. In 
our calculations it is satisfactory to look at the function values at x± « ±10, 
at least for the lowest eigenstates. The PT-symmetry of the potential plays 
now a key role. Since we are interested only in the real part of the spectrum 
we can assume the T'T-symmetry of the wave function. The equation for the 
eigenvalues, which originally readsQ 

is now, having in mind (|4ap and (|4b[). simplified into 

ReME,x+yME,x+)=0. (6) 

Note that the left hand side of ([5]) is a complex function of energy, whereas 
the left hand side of (O is real. This is clearly an advantage - the zeros of 
a real function can be determined by the bisection method. The described 
procedure works well for the case e = where comparison with the exact results 
is feasible; the exact result is recovered up to five- or six-digits precision, the 

•^This choice puts a limit on applicability of the method to e < 2. 

^The energy dependence of the solution is made explicit in the following. 
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Figure 4: Spectra of Bender Hamiltonian, contour with A = 1. The graph on 
the right hand side iUustrates the behaviour around e — 1 where real eigen- 
values emerge again. The algorithm was requested to find only lowest six real 
eigenvalues (five on the right graph), hence the lines corresponding to higher 
eigenvalues are discontinued where lower eigenvalue pair emerges. The chaot- 
ically distributed points in the upper part of the graph result from the error 
in the algorithm's implementation which occurs when the demanded number of 
real eigenvalues is not found. On both graphs, the vertical axis represents the 
energy while e is drawn horizontally. 
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Figure 5: A detail of the critical region between e = —0.45 and e — —0.65 for 
A= 1. 



precision can be enhanced with some loss of speed. Comparison of the non- 
tobogganic case with results obtained in [1] shows no significant differences. 
When £ — > 2 the computation becomes slower if the precision of computation 
has to be maintained since the solution's asymptotics becomes more vulnerable 
to numerical errors. 

The results show that the behaviour of eigenvalues depends on the winding 
number A. The results for A = 1 (see Fig. |3]) exhibit vast differences from 
the A = case. First, except the lowest eigenvalue the spectrum complexifies 
also at e > 0. As Fig. [3] suggests, the region of reality is broader for the 
low lying eigenvalues. It is entirely possible that there does not exist either 
left or right neighborhood of e = where the whole spectrum is purely real, 
in contradistinction to the non-tobogganic contour which yields real spectra 
for all e > 0. However it seems that there is a previously unattested region 
of real eigenvalues at £ < —0.4, probably not perturbatively accesible since 
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Figure 6: The A = 2 case. The region of reahty for each eigenvalue is narrower 
than for A = 1. On the right hand side the critical region is again depicted in a 
greater detail, its position is moved further to the right with respect to A = 1 
case. 



perturbative calculations usually break down in exceptional points. The lowest 
energy tends to infinity as£— >— linA = case; if A = 1 it joins the other 
real eigenvalues in the left region and eventually complexifies in an exceptional 
point near e = —0.61 (see Fig. [5|) - it is interesting to note that here it does not 
represent the ground state. The spectrum has to also be real in the vicinity of 
e G N since the singularity disappears there and the contours for different A are 
equivalent, consequently the real spectrum of A = case must be reproduced 
(see Fig. ED . 

The overall picture does not change significantly when A = 2. The overall 
pattern is similar to A = 1 (see Fig. [Sj . Possibly another region of real spectrum 
exists near e — —1, but the computation becomes lengthy and unreliable in those 
points, since near e = — 1 the solutions do not decrease enough rapidly; we are 
not confident in results obtained in this region by the above described method 
and leave this problem for future investigation. 



4 Summary 

The introduction of tobogganic contours into the Bender potentials produces 
another versions of PT-symmetric Hamiltonian. Though they are closely re- 
lated to the original non-tobogganic Hamiltonian, they are indeed different and 
exhibit qualitatively distinct behaviour with exceptional points standing be- 
tween intervals of reality including the points e = rt, n G N. In an alternative 
approach, one can change the variable to unbend the contour, this leaves e.g. 
the A = 1 Hamiltonian in the form 



11 _ ^ 4. ^%,6+3e 

9y4 y 



■fy'+^^. (7) 



after putting x — y'^ . Such transformations were discussed in [?]. They are in- 
teresting as a mehod that allows to simply transform the problem to an ordinary 
differential equation of one real variable. On the other hand the Schrodinger- 
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Figure 7: The first (left) and fifth (right) eigenvalue plotted against e for A = 
and A = 1 (the latter is the less straight dependence). 

like form of the Hamiltonian is lost, which makes the example less physcally 
appealing. 

If we are interested only in the transition between the harmonic oscillator and 
the negative quartic oscillator, the choice of A is clearly irrelevant. For integer 
e there is no singularity and distinctly winded contours must yield identical 
spectra. Therefore it can be said that any tobogganic Hamiltonian defines good 
continuation of the harmonic oscillator, and the A = special case is only 
"incidentally" privileged due to its real spectrum. 

It may be also noted that the dependence on A is non-perturbative in e. 
The equality of the linear approximation coefficients for distinct A is visible 
from Fig. [71 Up to the first order the energy is, independently of the contour 
selection, 
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